Effects on the eigenvalues of the quantum bouncer due to dissipation 
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^[ ABSTRACT 

a 

Effects on the spectra of the quantum bouncer due to dissipation are given when a 
linear or quadratic dissipation is taken into account. Classical constant of motions and 
^ '. Hamiltonians are deduced for these systems and their quantized eigenvalues are estimated 
■ through perturbation theory, we found some differences when we compare the eigenvalues 



of these two quantities. 
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I. Introduction. 

Dissipative systems has been one of the must subtle and difficult topics to deal with 
in classical [1] and quantum physics [2] . In general, to construct a consistent Lagrangian 
and Hamiltonian formulation for a given dissipative system can be a really challenge [3] . 
There are basically two approaches to study dissipative systems. The first one tries to 
bring about the dissipation as a results of averaging over all the coordinates of the batch 
system, where one considers the whole system as composed of two parts, our original con- 
servative system and the batch system which interacts with the conservative system and 
causes the dissipation (of energy) on it [4] . This approach has its own value and will not be 
followed or discussed here. The second approach considers that the bath system produces 
on our initially conservative system and average effect which is expressed as an additional 
external velocity depending force acting on the conservative system and transforming it 
into a dissipative system with this velocity depending force, the resulting classical dissipa- 
tive system contains then this phenomenological (or theoretical) velocity depending force. 
Then, the question arises over its consistent Lagrangian and Hamiltonian formalism and 
the consequences of its quantization. This approach, in addition, allows us to study and 
test the Hamiltonian approach for quantum mechanics and its consistence [5] and is the 
approach we will follow in this paper. A system which has bring our attention for dissipa- 
tion study through the above approach is the quantum bouncer. The quantum bouncer [6] 
is the quantization of the motion of a particle which is attracted by the constant gravity 
force, that is, close to surface of the earth. This particle hits a perfectly refiexing sur- 
face, producing the bouncing effect. This system with an additional dissipation force has 
particular importance because of its potential experimental realization. This dissipative 
system has been studied very little and using so far the first approach mentioned above 
[7]. 

In this paper we considers the second approach to dissipation systems and will assume that 
the external velocity depending force has linear and quadratic dependency with respect the 
velocity. This approach gives us the opportunity to check the nature of quantization via 
Hamiltonian or constant of motion associated to the system, that is, using the usual quan- 
tization of the linear generalized momentum or using the quantization of the velocity. This 
consideration is particularly interesting in dissipative systems since one can not always 
have a Hamiltonian in an explicit form as a function of the variables position and linear 
momentum [8], that is, the velocity "v" can not always be known explicitly in terms of the 
linear momentum "p" and position "x" of the particle through the relation p = dL/dv, 



2 



where L is the Lagrangian of the system. This paper is organized as follows: we present 
the classical study for the dissipative system considering the linear and quadratic velocity 
depending force. The constant of motion, the Lagrangian, and the Hamiltonian of the 
system are derived, and we give their expressions up to second order in the dissipation 
parameter. We presents the modification of the eigenvalues for the quantum bouncer, 
when this dissipation is taken into account, for the above approximated (weak dissipa- 
tion) constant of motion and Hamiltonian using quantum perturbation theory. Finally, 
conclusions and some discussions of our results are made. 

II. Classical linear dissipation. 

The motion of a particle of mass m under a constant gravitational force and a linear 
dissipative force is described by the equation 

= -mg - av , (1) 

where x is the position of the particle, g is the constant acceleration due to earth gravity, 
a is the parameter which characterizes the dissipation, and v = dx/dt is the velocity of the 
particle. A constant of motion of the autonomous system (1) is a function = K^ix, v) 
satisfying the equation [9] 

dKg a dKg 

V— {g+—v) = 0. (2) 

ox m V 

The solution of this equation such that limQ,_>o = mv'^ /2+mgx (the usual total energy 
of the non dissipative system) is given by 

TYh^ QV / TflQ \ ^ CXV 

Ka = — m — ln(l + — )+mgx. (3) 

a \ a / mg 

The Lagrangian associated to (1) can be obtained using the known expression [5] 

La = v — ^ dv , (4) 

bringing about the following Lagrangian 

m?qv ^ otv\ fmq\'^, av\ m^qv 

La = — — In lH +m — ln l + - mgx — . (5) 

a \ mg ) \ a / \ mg J a 

Therefore, the generalized linear momentum and Hamiltonian are given by 

m^g ( av 



p« = — In 1 + — (6) 
a \ mg / 
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and 



Ha = m[ — exp - 1 Pa + mgx . (7) 



At two orders in the dissipation parameter a, one has the constant of motion, the La- 
grangian, the generahzed hnear momentum, and Hamiltonian given as 

„ 1 2 a n A , 

K — —mv + mgx v H t^v , (8aj 

2 2g Amg-^ 

L = -mv — max v H t^v , (8o) 

p = mt; V H f , (8c) 

2g 3mg 

and 

2m 6mg^ 24m^g'^^ 

The constant of motion (3) or (8a) and the Hamiltonian (7) or (8d) bring about the 
damping bouncing effect on the spaces (x^v) and (x^p). The dissipative parameter a can 
be determined by measuring the velocity Vg at the refiexing surface {x = 0) and then 
measuring its maximum displacement x^ax {v = Q). Equaling the value of the constant 
of motion on both situations, one gets the expression 

- m — ln(l + — -) = mgxmax , (9) 

a \ a / mg 

from the parameter a can be gotten. 

III. Classical quadratic dissipation. 

In this case, the motion of the particle is described by the equation 

tJC' II . 

TO-^^ = —mg — ^v\v\ , (10) 

where 7 represents a dissipation constant which, of course, is different from the previous 
case. Proceeding in the same way as we did for the linear case, the constant of motion, 
Lagrangian, generalized linear momentum, and Hamiltonian are given by 
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= mvexp I ± I , (lie) 



m 



and 



where the upper sign corresponds to the case v > 0, and the lower sign corresponds to 
the case v <0. These equation where aheady given in reference [10]. The damping effect 
of the bouncing particle in the space (x, v) can be traced in the following way: starting 
with the initial condition Xo = and Vo > 0, for example, the constant of motion is 
determined, = mv'^/2. Then, the maximum distance Xmax {v = 0) is calculated from 
the expression = (m'^g/2'y){exp{2'jXmax/^) — i) which helps to calculate the constant 
K-^ K- = —{m'^g/2'y){eyi-p{—2'^Xjnax/^) — l)- This K- is used now to calculated the 
velocity at the return point {x = Q)^ v\ = —\j2K-/'m. Due to perfectly reflexing surface, 
the velocity of the bouncing particle for the next cycle is vi = —v^ {vi < fo), and the 
above cycle is reproduced again, and so on. Starting with the same initial conditions, the 
trajectories in this space are one below the other at any time, as the damping factor is 
grater. The damping effect in the space (x, p) through the Hamiltonian approach can be 
analyzed similarly. However, the trajectories starting with the same initial conditions on 
this space are not below the other all the time, as the damping factor is greater. This 
strange effect is due to change in sign in (lid) with respect to (11a), produced by the 
position and velocity dependence of the expression (11c). 

To determine the constant 7 through the constant of motion, one can start with the 
initial conditions {xg = 0, > 0) and can determined the constant of motion = 
mVg/2. Then, one can measure the maximum displacement Xmax {v = 0) and to solve 7 
from the equation 

''„ = '4l(e.J?2f^)-l). (12) 



2 " 27 V V ™ 

Up to second order in the dissipation parameter, one has from (11a) to (lid) the constant 
of motion, the Lagrangian, the generalized linear momentum, and the Hamiltonian given 
by ^ 

K± = -mv^ + mgx ± j[v^x + gx^] + ^\v^x^ jm + 2^a:^/3m] , (13a) 

L± = ^mv'^ — mgx ± 'y[v'^x — gx^] + ^\v^x^ jm — 2gx^ /2im] , (13^) 
= mv ± 7[2t;a:] + ^'^[2vx'^ /m] , (13c) 
iy± = |— + mgx =F ^[p^x/m^ - gx'^] + 7^[p^a;^/m^ + 2gx^ /im] . (13d) 
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IV. Quantization of the constant of motion. 

Eqs. (8a) and (13a) can be written as 



K{x, v) = Ko{x, v) + V{x, v) , 



(14) 



where Ko is the constant of motion without dissipation 

Ko{x, v) = -mv^ + mgx , 
and V takes into account the dissipation factors 



(15a) 



V{x, v) 



^7 [v'^x + gx'^] + 7^ 



IZV , 2gx'' 



(linear case) 
(quadratic case) 



(156) 



The quantization of (14) can be carried out through the associated Schrodinger's equation 
of this constant of motion 

(16) 



where ^ = ^{x, t) is the wave function, h is the Plank constant divided by 27r, K = Ko+V 
is a Hermitian operator associated to (17), and v is the velocity operator defined as 



^ ih d 
m dx 



(17) 



Since Eq. (16) represents an stationary problem, the usual proposition 
'if{x,t) = exp{—iE^t/h) 'il){x) transforms (16) to an eigenvalue problem 



(18) 



Considering now the operator F as a perturbation of the constant of motion Ko-, one can 
calculate an approximated solution to the problem (18) through perturbation theory . 
The solution of the eigenvalue problem 



(19) 



is well known [6] , with il^n^ being the eigenfunction given by 



Ai{z - Zn) 
\M'{-Zn)\ 
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(20) 



where Ai and Ai' are the Airy function and its differentiation, and Zn is its nth-zero 
{Ai{—Zn) = 0) which ocurres for negative argument only, z is the normalized variable 
z = x/lg with Ig = (h^ /2m'^g)^^^ , and z^ is related with the eigenvalue through the 
expression 



Zr}. 



mglg 



(21) 



Up to second order in perturbation theory, the eigenvalues of (18) are given (in Dirac 
notation [11]) as 

{n\V\k)\'^ 



E:^ = E^'^ + {n\V\n) + J2 



(0) 



(22) 



where {z\n) 



Using the Hermitian operators v'^x = (v'^x + vxv + xv"^)/?) and 



v^x 



V X + vx V + X V + XV X + xvxv + vxvx)/6 for the associated expressions on 



(25b), and using the relations {n\x^\k) = lg{n\z''^\k) and {n\d'^ / dx^ \k) 
for any integer s, one has (see appendix for a list of matrix elements) 



E^ = Ei'^ + l 



a 



8^,3^|l/2 + m<7/,/(£;r-£;,f )| 



k^n 



(0) 



T 7 



15 +^ 



E 



1 564 
"5 + TDS" 



2gl 



m 



where a„fc are real number given by 

1 12 - 2zk{Zn - Zkf + {Zn - Zk)'^\^ 



(^nk 



{Zk - Zn)" 



l-'{n\d'/dz'\k) 

(linear) 

, (quadratic) 
(23a) 

{23b) 



Note that for the linear dissipation case, there is not real contribution a first approxi- 
mation, and for the quadratic dissipative case, the first order contribution depends on 
whether the particle is moving (-) or down (-). Within a full cycle, this first order is can- 
celled out and the second order contribution remains. Of course, for the approximation 
(3a) to be valid, one must have that the second term in this expression must be much less 
than En"* which restricts the value of the dissipative parameter. 
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V. Quantization of the Hamiltonian. 

Eqs. (8d) and (13d) can be written as 



H{x,p) = Ho{x,p) + W{x,p) 



(24) 



where Hq is the Hamiltonian without dissipation, 



p 

Ho{x,p) = \-mgx , 

2m 



(25a) 



and W has the dissipation terms, 



W{x,p) = 



( 

-r- \P^^ 

m 



+ a' 



p 



5„2 



24m°g 



2 2 r> 

2fP X 2gx'- 



m 



(hnear) 
(quadratic) 



(256) 



It is necessary to mention that the quantization of some systems for quadratic dissipation 
has been solved by different authors [10,12] but perfectly reflexing wall potential. 



V{x) = 



f oo for a: < 



mgx for a: > 



(26) 



Moreover, the solution given in reference [10] is singular for the dissipation parameter 
equal to zero. Therefore, we think it worths to make the analysis of the quantization for 
small orders in the parameter 7. For the usual Shrodinger quantization approach, one has 
the stationary equation 

(27) 



ih— = H{x,p)^ 



where H is the Hamiltonian operator associated to (24), and p is the usual linear mo- 
mentum operator p = —ihd/dx. Eq. (27) is transformed to an eigenvalue problem, 
H'ip{x) = E^ip{x), through the proposition "^(xjt) = exp (^—iE^t/h) ijj(x). Since the 
Hamiltonian H is given by H = Hg + where the solution of the equation 



(28) 



is well known (as before), being il^n^ and En^ given by (20) and (21), perturbation theory 
can be used to determine the approximated values of the eigenvalues E!^^ (similarly as 
done with expression (22)). Using the Hermitian operators p'^x = {p^x -\-pxp-\- xp^)/S 



8 



and p^x^ = {p^x^ + px^v + x'^p^ + a;p^a; + xpxp + pxpx)/6 for the associated expressions 
on (25b), one gets 



a 



-4°' 



(linear) 



2^2 



1 ^ 564\ Ml^/i 2,4 



k^n 



, (quadratic) 
(29) 

where a^k is given by (23b). As one can see from (23) and (29), there is a difference 
with the eigenvalues associated to the constant of motion and Hamiltonian quantizations. 
Their relative differences, 5En = {E^ — E^)/En \ is given by 



a 



= < 



1 IgZ. 



k^n 



|l/2 + mg/,/(i^f -EW)| 

^(0)_^(0) 



(linear) 



(quadratic) 



(30) 



VI. Conclusion. 

We have presented the classical and quantum analysis of a particle attracted by constant 
gravity and velocity depending dissipative forces which bounce on a perfectly reflexing 
surface. We have considered linear and quadratic velocity depending dissipative cases and 
have deduced their constant of motion and Hamiltonians. Expression (3) and (11a) gives 
us the expected damping behavior of the particle on the space {x,v), but the expression 
(7) and (lid) show us the damping with unexpected behavior in the space {x,p) (two 
trajectories on this space, for dissipative parameter one bigger than other, do not follow 
one under the other all the time). For the quantum case, we have analyzed the eigenvalues 
of the assigned operators pertubedly and up to second order on the dissipative parameters. 
Relation (30) tells us that that there is a difference whether the constant of motion or the 
Hamiltonian is quantized, and it suggests that one could see this difference experimentally. 
In this way, one could see whether nature prefers to follow constant of motions rather 
than Hamiltonians for dissipative systems. Finally, one must observe that for the full 
linear case (7), it is possible to solve exactly the Shrodinger equation in the momentum 
representation, and this will be analyzed on a future paper. 
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APPENDIX 



We show here a hst of some matrix elements from reference [6] and some other calculated 
from the same reference (an correction of a sign has been made to some matrix elements). 
Given the functions (20) and n 7^ A;, one has 

{n\k) = 5nk {Ai) 



{n\z\n) = -Zn {n\z\k) = {A2) 

^ [^n ^k) 

/ I 3| \ 3 48 3 3 24(^^ + ^fe)(-l)"+fe+^ 

{n\z \n) = - + —z^ {n\z \k) = ^^---^^ (^4) 

{n — n)=Q {n — k) = ^ 

dz dz Zn — Zk 

/ , , V 1 , id\,^ 2(-l)"+^ 

{n\i,\n) = I {n\l,\k) =(l + (-1)^+^ (.1.) 



dz^ 2 dz^ V 2 Zk — z> 



1 

/ I " I \ 2 



{n\ — \k) = ~2(^fc - ^n) + 24 - 2zkizk - Znf ^ 



'dz^ 



{Zk - ZnY 
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